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The complex Langevin method is a leading candidate for solving the sign problem occurring in 
various physical situations, notably QCD at finite chemical potential. Its most vexing problem is 
'convergence to the wrong limit', where the simulation gives a well defined, but incorrect, result. 
Here, we first outline a formal justification of the method and identify points at which it might 
fail. From these we derive a condition that must be satisfied in order for correct results to be 
obtained. We then apply these ideas to the three-dimensional SU(3) spin model at finite chemical 
potential and show strong indications that complex Langevin dynamics yields correct results in 
this theory. 
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1. Introduction 

Many systems of physical relevance have a complex action, QCD with a finite chemical poten- 
tial is a particularly interesting example. Such theories are difficult to study numerically because 
the weight is complex and therefore standard Markov chain Monte Carlo techniques based on a 
probability interpretation cannot be applied; this is often called the sign problem [|l|]. The complex 
Langevin equation (CLE) enjoyed a certain limited popularity after being proposed in the early 
1980s ^ ^] because it can, in principle, avoid the sign problem since it does not rely on impor- 
tance sampling. In some cases, complex Langevin simulations have been found to be numerically 
unstable due to runaway solutions but this can be cured by use of an adaptive stepsize [^. A more 
serious and vexing issue is convergence to the wrong limit, where the simulation gives well defined 
results, but when compared with known values are found to be incorrect. It is this problem that we 
shall focus on here. 

In Section |2| we briefly cover a formal argument for the correctness of the CLE and identify 
points at which it might fail. By studying the long time evolution of observables with respect 
to real and complex measures, we derive a criterion which must be satisfied in order for correct 
results to be obtained [^, ^]. In Section ^ these ideas are applied to the SU(3) spin model at finite 
chemical potential ||7|, |8|, ^. A comparison with imaginary chemical potential shows that the CLE 
works correctly in the region jj.^ ~ 0. An analysis of the criterion with a larger chemical potential 
supports the claim that in this case complex Langevin dynamics can be trusted [p^]. 



2. Criteria for correctness 



The central object of interest is the expectation value of a particular observable, given by 

J 0{x)e-^^''Ux 



(O) 



(2.1) 



where for notational simplicity we use a single real degree of freedom, x. The action S{x) is 
complex, preventing a probability interpretation of the measure and ruhng out methods based on 
importance sampling. The complex Langevin equation is 



dx _ 



/NrTIr, ^=K, + ^Nir]i, (2.2) 

where the real variable is analytically continued as x — z = x + /y. The drift terms are given by 



dS{x) 



dx 



dS{x) 



Ky = \m ^ 
dx 



(2.3) 



and the two noise terms tjr, T]i are independent Gaussian random numbers with variance 2 and nor- 
malisation A^i > and N^—Ni = 1. A numerical simulation can then be implemented by integrating 
these equations to large times f — oo. 

The resulting dynamics is described by a dual Fokker-Planck equation for the evolution of the 
probability density P{x,y;t), 

^P{x,y;t)=L^P{x,y;t), (2.4) 
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with the operator 

= V, [A^rV, - K,] + VylNiVy - Ky] . (2.5) 

To understand the time evolution of the real density P{x,y;t) one must also examine the evo- 
lution of the complex density p{x;t), determined by 

^^p{x;t)=Llp{x;t). (2.6) 

Here, the complex Fokker-Planck operator Lq is 

4 = V.[V, + V,5(x)]. (2.7) 

This equation has p(x;oo) oc exp[— 5(x)] as a stationary solution, which is expected to be unique. 
Numerical studies (where feasible) of Eq. (^!^ ) confirm this to be true; in fact, convergence to this 
distribution seems exponentially fast. 

Expectation values with respect to the two densities can now be defined as 



in\ _ IO{x + iy)P{x,y;t)dxdy _ J 0{x)p{x;t)dx 

^^^^W - jP{x,y;t)dxdy ' ^^^''^ " fp{x;t)dx ' ^^"^^ 



The result that one would like to show is 

{0)pit) = {0)pir), (2.9) 

if the initial conditions (C?)p(o) = (O)p(o) match, which is assured provided 

P{x,y;0)=pix;0)5iy). (2.10) 

One expects the dependence on the initial conditions to vanish in the limit f — )• oo by ergodicity. 

To establish a connection between the expectation values with respect to P and p, one moves 
the time evolution from the densities to the observables. Since we are only interested in functions 
of z = x + iy (holomorphic functions), we may act with the Langevin operator 

whose action on holomorphic functions agrees with that of L. 

We now use L to evolve observables according to the equation 

^0{z;t)=LO{z;t), (2.11) 

which is formally solved by 

=exp[fL]0(z). 

Due to the fact that L and L agree on holomorphic functions, the tilde may be dropped. 
To examine the evolution we define the function 

F(?,t) = J P{x,y;t-T)0{x + iy;T)dxdy, (2.12) 
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and observe that F(t, t) interpolates between the two expectation values: 

F{t,0) = {O)p^,), F{t,t) = {0)p(^,y (2.13) 
The first can be seen easily, while the second makes use of the initial conditions and 
F{t,t) = J P{x,y;0) [e'^O) {x + iy;0)dxdy 

= j p(x;0) [e'^O] {x;0)dx = J O(x;0) (/^"p) ix;0)dx = {0)p^,y (2.14) 

where it is only necessary to assume that integration by parts in x does not produce any boundary 
terms. 



The desired result (2.9) follows if F{t,T) is independent of T. To check this, we need the T 
derivative to vanish. 



J {L^P{x,y\t-z))0{x + iy\z)dxdy + j P{x,y;t - z)LO{x + iy\z)dxdy. (2.15) 



Integration by parts, if applicable without boundary terms at infinity, then shows that the two terms 
cancel and therefore F{t, x) is independent of T, irrespective of N\. 

This is therefore a point at which the formal argument might fail: if the decay of the product 

P{x,y,t -T)0{x + iy\T:) 

and its derivatives is insufficient for integration by parts without boundary terms. 

In Ref. 1^ a study of the U(l) one-Unk model found that the T-derivative is largest at T = 0. 
This motivates the superficially weaker condition 



lim —F(t,T:) 
dt 



= 0. (2.16) 

T=0 



This modification to the condition is still sufficient for correctness if it holds for a sufficiently large 



set of suitably chosen observables. Taking the limit ? — oo in Eq. (2.15) causes the first contribu- 
tion to vanish because of the equilibrium condition, P{x,y;oo) = 0. Therefore, the criterion for 
correctness reduces to 

Eo = J Pix,y;oo)LO{x + iy,0)dxdy = {LO) = 0. (2.17) 
This is fairly simple to check for a given observable, but it is in fact a strong statement since it 



must hold /or all observables. Therefore, Eq. (2.17) really represents an infinite tower of identities 



which must all be satisfied. In practice it can be checked for a small number of observables, still 
yielding a necessary criterion [^]. 

3. SU(3) spin model 



Motivated by recent work [|9p (for related models see e.g. Ref. []1 1[]), we re-examine the three- 
dimensional SU(3) spin model at finite chemical potential, for which promising results with com- 
plex Langevin dynamics have been obtained in earher studies [^, ||]. The action is formed from 
three contributions 

S = Sb + Sf + Sh, (3.1) 
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which read 



Sb = -j8£ (Tr[/,Tr[/7^l +Tr[/7iTr[/,+,) , 

x.v 

Sp = {e^TrU:, + e-^TrU^^) , 

X 

5H = -£log 



Sin - — sin'' — — . sin-' ' 



(3.2) 
(3.3) 

(3.4) 



The final component, Sh, originates from the Haar measure introduced by diagonalising U in terms 
of angles, 

The action is complex due to the "heavy fermion" contribution, 5'p(/x) = 5'f(— /^*). For small 
h the theory has a confined and deconfined phase separated by a first order transition line at small 
jj., turning into a crossover at larger /i (see Figure [T], left) [^. 




Figure 1: Left: sketch of the phase diagram at small h. Right: (Tr {U + U ' ) /2) as a function of /3 at /i = 0, 
h = 0.02 on a 10^ lattice. 

Figure |l| (right) shows the transition as a function of j8 at = 0, using the observable (Tr ([/ + 
U^'^)/2). Note that at /i = 0, (Tr?7) = (Tr but at /^^ > they differ. Since observables that 
are even in pL should be continuous across /x^ = 0, we also perform simulations with an imaginary 
chemical potential, for which the action is real and simulations using standard techniques can be 
used (we use real Langevin dynamics). A comparison with complex Langevin data should therefore 
show even observables to be continuous across jx^ = 0. In Figure^ the observable (Tr(?7 + U^'^)/2) 
is plotted. The data from complex Langevin dynamics with /i^ > is consistent with those from 
real Langevin dynamics with /x^ < in both phases, including the critical region, suggesting that 
complex Langevin dynamics is working correctly here. This is in contrast to the case of the XY 
model recently studied using complex Langevin dynamics, where correct results were obtained in 



only part of the phase diagram [12|. 

At larger values of /i^, one can no longer rely on analytical continuation to justify the results. 
Instead, we assess them using the criteria developed above and compute {LO) for O = TrU . The 
outcome is shown in Figure ^ for real chemical potential (left) and also in the phase-quenched 
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Figure 2: (Tr(f/ + f/-i)/2) as a function of for various j3 values. For /z^ < (imaginary /i) the action 
is real and real Langevin dynamics is used; when /i^ > complex Langevin is necessary. 
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(a) Complex Langevin dynamics 



(b) Phase-quenched using real Langevin dynamics 



Figure 3: Langevin stepsize dependence of the correctness criterion for Trt/ with complex Langevin dy- 
namics (left) and in the phase-quenched approximation (right) using real Langevin dynamics. The points for 
\i = 3.5 are shifted horizontally for clarity ()3 = 0.12, h — 0.02 and lattice volume 4^). 



theory (right), using /i = 3 and 3.5. In the phase-quenched theory, the action is real and we use real 
Langevin dynamics. The figures indicate that at nonzero stepsize the criterion is not satisfied, both 
for real and complex Langevin dynamics, but that in the limit e — )^ it is satisfied. 

We find therefore that in this model the test is passed successfully. Moreover, finite stepsize 
corrections can be quantified by the deviation of {LO) from zero. We are currently extending 
this analysis and have implemented a higher-order algorithm to eliminate the effects of finite step- 
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size [0]. 

4. Conclusions 

Complex Langevin dynamics can in principle be applied where the sign problem prevents the 
use of importance sampling. An analysis of the long-time evolution of real and complex measures 
using the Fokker-Planck equation shows that the correct stationary solution exists. However, if 
the decay of the distribution P{x,y;t) is insufficient to allow integration by parts without boundary 
terms, convergence to the wrong limit can occur. In practice, this can be diagnosed by testing the 
condition that (LO) = for a suitably large set of observables in the limit of vanishing stepsize. In 
the case that the criterion is not satisfied, complex Langevin dynamics fails. 

An analysis of the SU(3) spin model shows that the criterion for correctness is satisfied, justi- 
fying the claim that complex Langevin dynamics works correctly with this model. This is corrob- 
orated by analytic continuation from imaginary chemical potential in the region of small chemical 
potential |jU^| < 1, which shows that data for observables even in /i are continuous between results 
from a real action when /i^ < and a complex action with > 0. A further and more detailed 



study of the criterion and stepsize dependence is in progress [|10|. 
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